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1. Prologue 

3D Navier-Stokes equations (NSE) - describing a flow of 3D incompressible vis- 
cous fluid - read 

Ut + (u ■ V)u = — Vp + Am, 
supplemented with the incompressibility condition div u = 0, where u is the ve- 
locity of the fluid and p is the pressure (here, the viscosity is set to 1); taking the 
curl yields the vorticity formulation, 

uj t + (u ■ V)w — (uj ■ V)u + Alo, 

where cj = curl u is the vorticity of the fluid. 

It is well known that both globally [ Le34 . H51 [ and unif ormly-locally (with suit- 
able spatial decay at infinity) I L-R02I finite energy data generate global-in-time 
weak (distributional) solutions to the 3D Navier-Stokes equations (NSE), satisfy- 
ing global and local energy inequality, respectively. Despite much effort - since 
the pioneering work of Leray [Le34[ in 1930's - the question of whether weak so- 
lutions may exhibit finite-time singularities remains an open problem. It is known 
that the set of all possible singularities is small - the one-dimensional (parabolic) 
Hausdorff measure of the singular set in VL x (0, T) is zero for any T > [CKN82|; 
here, Q is a global spatial domain. 

There are various regularity criteria preventing finite-time formation of singu- 
larities, mainly expressed either as a local or a global condition on a weak solution 
over a spatiotemporal domain, or as a condition on a regular solution approaching 
a potential singular time T* . The conditions are given as boundedness of a suitable 
spatiotemporal norm, the common trait being that the norm is scaling-invariant 
(critical) with respect to the natural scaling in the Navier-Stokes model. In con- 
trast, the a priori bounded quantities are all subcritical; moreover, there is a scaling 
gap between a regularity criterion in view and the corresponding a priori bounded 
quantity. As an illustration, here are two classical examples - for the velocity 
formulation - in L p and BMO spaces. The regularity criteria are boundedness 
in L^°L^ [ISS03J and L\BMO x [KTOOj, and the corresponding a priori bounded 
quantities are Lfl? x HLe34||H51| and L\Lf HFGT81I , respectively. These are man- 
ifestations of supercriticality of the Navier-Stokes problem. 

A rigorous study of geometric depletion of the nonlinearity in the 3D NSE (as well 
as in the 3D Euler equations) was pioneered by Constantin, and it was based on 
a singular integral representation for a - the stretching factor in the evolution of 
the vorticity magnitude \uj\ depleted by local coherence of the vorticity direction - "the 
story of alpha and omega" [Co94[. This is fundamental as there is ample evidence 
- both numerica l IIAKKG87I [IWSK93J I58T1 1 VM941 gTOgl) and theoretical IICPS951 
IGGH97[ IGFD991 |Oh09 1 - that regions of intense vorticity tend to self-organize in 
coherent vortex structures, most notably, quasi one-dimensional vortex filaments, 
displaying strong local coherence of the vorticity direction. 

The mechanism of the geometric depletion of the nonlinearity was subsequently 
exploited in [ CoFe93 1 to show that as long as the regions of intense vorticity ex- 
hibit local Lipschitz-coherence of the vorticity direction, no finite-time blow up 
can occur, and later in [ daVeigaBe02 [ where the Lipschitz-coherence was replaced 
by the ^-Holder coherence. A spatiotemporal localization of the ^-Holder coher- 
ence regularity criterion was performed in [GrZh06 Gr09|, and independently in 
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BChKaLe07l . The aforementioned regularity criteria are all pointwise coherence 
conditions; hence, necessarily supercritical with respect to the NSE scaling. A lo- 
cal, scaling-invariant (critical) criterion over a parabolic cylinder below a potential 
singular point (xo, to), 



\ui(x,t)\ pi 2 „(i, t) dx dt < 



- 2R 

'*o-(2i?) 2 JB(x Q ,2R) 2 ' 

where 

smcp(r](x, t),r)(y,t) 



p~/^ r (x, t) = SUp 



yeB(x,r),y^x \x ~ Vl 1 

is a 7-Holder measure of coherence of the vorticity direction rj at the point (x, t), 
was presented in IGrG ulCPTl . On the other hand, a corresponding (subcritical) a 
priori bound had been previously obtained in |Co90L 

\u)(x,t)\\\7ri(x,t)\ 2 dxdt < -v~ 2 I \u (x)\ 2 dx 

2 .Id3 



U 



where v is the viscosity. 

A different geometric approach to the study of possible singularity formation 
in 2D and 3D incompressible flows was developed by Cordoba and Fefferman 
nCF01-l[|CF01-21ICF 02l; in particular, non-existence of 'tube collapse singularities' 
in 3D incompressible inviscid flows was shown in |CF01-2], and non-existence of 
a more general class of 'squirt singularities' in incompressible flows - including 
the flows described by the 3D NSE - was presented in |CFD04). 

The purpose of this Note is to present a mathematical evidence - in a sta- 
tistically significant sense - of a geometric scenario leading to criticality of the 
Navier-Stokes problem. More precisely, utilizing the ensemble averaging process in- 
troduced in our recent study of turbulent cascades in physical scales of 3D incom- 
pressible fluid flows HDaGrll-11 IDaGr 11-21 IDaGrl2-ll IDaGr 12-21 , we show that 
the ensemble-averaged vortex stretching term is positive across a range of scales 
extending from a square root of a Kraichnan-type micro-scale to the macro-scale. 
Combining this with the a priori bound on the decrease of the distribution func- 
tion of the vorticity obtained by Constantin in [Co90[ - as well as with the general 
mechanism of creation and dynamics of vortex filaments in turbulent flows (cf . 
[CPS95Q - indicates a geometric scenario in which the region of intense vorticity 
(defined as the region in which the vorticity magnitude - near a possible singular 
time - exceeds a fraction of the L°°-norm) comprises of macro-scale-long vortex 

filaments with the diameters of the cross-sections scaling like — . This is 

ll"(t)IIS> 

exactly the scale of local one-dimensional sparseness of the region of intense vorticity 
needed to prevent a formation of a finite-time singularity BGrl2l . 

2. Geometric measure-type regularity criterion 

In this section, we briefly recall a geometric measure-type regularity criterion for 
solutions to the 3D NSE obtained recently by one of the authors; for details, see 
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Definition 2.1. Let xq be a point in R 3 , r > 0, S an open subset o/R 3 and 8 in (0, 1). 

The set S is linearly 5-sparse around xq at scale r in weak sense if there exists a 
unit vector d in S 2 such that 

\S n (xq — rd, xq + rd)\ 



2r 



< 8. 



For M > 0, denote by f2 t (M) the vorticity super-level set at time t; more pre- 
cisely, 

n t (M) = {x e R 3 : \ui(x,t)\ > M}. 
The vorticity version of the local one-dimensional (linear) sparseness regularity 
criterion is as follows. 

Theorem 2.1. Suppose that a solution u is regular on an interval (0, T*). 

Fix 8 in (0, 1), and let h = h(5) = -| arcsin j^jp and a = a{8) > Assume that 
there exists e > such that for any t in (T* — e, T*), either 

(i) t + ,„,. 7 ... — > T* (do is an absolute constant appearing in the local-in-time 
analytic smoothing in L°°; cf. IGrl2|), or 

(ii) there exists s — s(t) in t + 4d a^ t ^ , t + ^jj^( t )|| such that for any spatial 
point xq, there exists a scale r, < r < j-, with the property that the super-level 

2dg|k(t)||J 

set fl s (M) is linearly 8-sparse around xq at scale r in weak sense; here, M = M(8) = 



Then, there exists 7 > such that uj is in L°°{JT* - e, T* + -f);L°°j, i.e., T* is not 
a singular time. 

The proof is based on a very intimate interplay between the diffusion in the 
model - represented by the local-in-time analytic smoothing in L°° - and the geo- 
metric properties of the harmonic measure (via the harmonic measure majorization 
principle). 

The analyticity estimate on solutions needed is a vorticity version of the esti- 
mate given in IGulOl ; this was based on a general method for estimating uniform 
radius of spatial analyticity in L p -spaces introduced in |GrKu98|, which was in 
turn inspired by the (analytic) Gevrey-class method presented in | FT89 ] (see also 
IIFeTi98l ). 

The key geometric harmonic measure estimate used in the proof is a generaliza- 
tion of the classical Beurling's problem |Beu33], conjectured in |Seg88| and solved 
by Solynin in | Sol99 1 (a symmetric version of the problem was previously resolved 
in |EssHa89|); for more details, see |Grl2). 

Remark 2.1. A rudimentary version of Theorem 12.11 was previously obtained in 
BGrOH . The condition needed in |Gr01] is a much stronger condition; essentially, 
a requirement of a local existence of a sparse coordinate projection. In contrast, all 
that is needed here is a local sparseness of an one-dimensional trace in a very weak 
sense. 
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3. The region of intense vorticity 

There is strong numerical evidence that the regions of high vorticity organize 
in coherent vortex structures 155TI IAKKG871 [SJQ91J [TWSR93J |VM94| and in par- 
ticular, in elongated vortex filaments (tubes). In addition, an in-depth analysis of 
creation and dynamics of vortex tubes in 3D turbulent incompressible flows was 
presented in ICPS95I (see also IGGH97llGFD99UOh09» . 

Consider a flow near the first (possible) singular time T* , and define the region 
of intense vorticity at time t < T* to be the region in which the vorticity magnitude 
exceeds a fraction of Ho^t)!^; keeping the notation from the previous section, this 

corresponds to the set fi t ^ — ||o»(t) ||oo for some c% > 1. 

Denote a suitable macro-scale associated with the flow by i? • The picture painted 
by numerical simulations indicates that the region of intense vorticity comprises - 
in a statistically significant sense - of vortex filaments with the lengths comparable 
to R Q . 

Let us for a moment accept this as a probable geometric blow up scenario. The 
length scale associated with the diameters of the cross-sections can then be esti- 
mated indirectly, by estimating the rate of the decrease of the total volume of the 

region of intense vorticity fi t ( — ||w(i)|| 

Taking the initial vorticity to be a finite Radon measure, Constantin showed 
[Co90| that the L 1 -norm of the vorticity is a priori bounded over any finite time- 
interval; a desired estimate on the total volume of the region of intense vorticity 
follows simply from the Tchebyshev inequality, 

voi(a(^iM*)iu))< ]j J c (c 3 >i). 

This implies the decrease of the diameters of the cross-section of at least 



IM*)IJ~ 

(03 > 1), which is exactly the scale of local one-dimensional sparseness of the region 
of intense vorticity needed to prevent the formation of singularities presented in 
Theorem 12. II In other words, the Navier-Stokes problem in this scenario becomes 
critical. 

A key step in justifying this scenario is providing a mathematical evidence of per- 
sistence - in a statistically significant sense - of the i?crlong vortex filaments (at 
this point, the evidence is purely numerical). A term responsible for the creation 
of vortex filaments is the vortex-stretching term, 

(to ■ V)u • to = Slo ■ to, 

where S is the strain matrix. One way to identify the range of (longitudinal) scales 
at which the dynamics of creation and persistence of vortex filaments takes place 
is to identify the range of scales ofpositivity of 5a; • to. In the following section, we 
will show that the range of positivity of Sco ■ to - in a statistically significant sense 
- extends from a power of a Kraichnan-type micro-scale to the macro-scale Rq . It 
is worth pointing out that the argument is dynamic - based on ensemble averaging 
local dynamics described by the full 3D Navier-Stokes system. 
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4. A DYNAMIC ESTIMATE ON THE VORTEX-STRETCHING TERM ACROSS A RANGE 

OF SCALES 

We begin by recalling the concept of ensemble averaging with respect to (K\ , K 2 )- 
covers at scale R, introduced in our work on existence and locality of turbulent cas- 
cades in physical scales of 3D incompressible flows 1 Da Grl l-lllDaGr 1 l-2UDaGr 12-11 
IDaGrl2-2H (for more details, see, e.g., HDaGrl2-2ll ). 

Let Rq > 0, and assume (for convenience) that the macro-scale domain of inter- 
est is the ball B(0, R ), B(0, 2R ) contained in the global spatial domain Q. Con- 
sider a locally integrable physical density of interest /, and let < R < Ro; the 
time interval of interest is (0, T). 

In what follows, we utilize refined spatiotemporal cut-off functions <p = 4>x ,R.t = 
%j} n, where v = r) T {t) € C°°(0,T) and ip = ip Xo , R (x) e V(B{x ,2R)) satisfying' 

(4.1) < r? < 1, r/ = 0on(0,T/3), r, = 1 on (2T/3, T), < ^ 
and 

(4.2) 0<*<1, t-lcBfr.*,, M<§, ^<%. 

for some | < pi, p 2 < 1. In particular, O = ^oV where ipo is the spatial cut-off (as 
above) corresponding to x = and R = R . 

For xq near the boundary of the macro-scale domain, 5(0, Rq), we assume ad- 
ditional conditions, 

(4.3) < V < "00 

and, if B(x ,R) <£ B{0,R o ), then tfj e V{B{0,2R )) with i/> = 1 on B(x ,R) n 
B(0, Rq) satisfying, in addition to l|4.2|l , the following: 



(4.4) 
and 



if) = ipo on the part of the cone centered at zero and passing through 
5(0, Ro) n B(x , R) between 5(0, Rq) and 5(0, 2R a ) 



tp = on B(0, R ) \ B(xo,2R) and outside the part of the cone 
(4.5) centered at zero and passing through 5(0, Ro) n B(xq, 2R) 
between 5(0, Rq) and 5(0, 2R ). 

A physical scale R is realized via suitable ensemble averaging of the localized 
quantities with respect to '{K\ , K 2 ) -covers' at scale R. 

Let Ki and K 2 be two positive integers, and < R < R . A cover {B(xi, i?)}™ =1 
of the macro-scale domain B(0, R ) is a (K\, K 2 ) -cover at scale R if 

^° ^\ 3 ^ s v ( r° 

-Rj ^ n ^ Kl {R 

and any point x in B(0, Ro) is covered by at most K 2 balls B(xi, 2R). The param- 
eters K\ and K 2 represent the maximal global and local multiplicities, respectively. 
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Considering the time-averaged, per unit mass - spatially localized to the cover 
elements B(xi,R) - local quantities f Xi ,R, 

fxt,R=7f; 7^/ f{x,t)(/) 6 xRT (x,t)dxdt 

1 JO H Jb(x,,2R) 

(for some < 6 < 1), the ensemble average {F)r is defined as 

1 n 

(F) R = - } f Xi ,R ■ 

n 

i=l 

The ensemble averages (with the fixed multiplicities K\ and K 2 ) act as a 'de- 
tector ' of significant sign-fluctuations of the density in view. More precisely if the 
density exhibits significant sign-fluctuations on the scales comparable or greater 
than R, the ensemble averages at scale R - with respect to all admissible (Ki, -Re- 
covers - will respond by exhibiting a wide range of values, from positive through 
zero to negative. This can be seen by rearranging the cover elements to emphasize 
the positive and the negative parts of the density, respectively. The larger the mul- 
tiplicities, the finer the detection. In contrast, if the ensemble averages at scale R - 
with respect to all admissible (K 2 , K 2 )-covers (again, with the fixed multiplicities) 
- are nearly independent on the particular choice of the cover, and say positive, 
this indicates that the density is essentially positive on the scales comparable or 
greater than R. 

As expected, for a non-negative density /, all the averages are comparable to 
each other throughout the full range of scales R, < R < Rq; in particular, they 
are all comparable to the simple average over the integral domain. More precisely, 

(4.6) -^F < (F) R < K*F 

for all < R < Ro, where 

and = K*(K U K 2 ) > 1. 

Consider now a global-in-time weak solution u (say, a global-in-time 'local 
Leray solution' on R 3 x (0, 00) in the sense of IL-R02J), and let T be the first (pos- 
sible) singular time. 

A spatiotemporal localization of the evolution of the enstrophy was presented 
in MGrZh06l IUr09|. Considering a (Ki, K 2 ) -cover {B(x l , at scale R, the 

following expression for the time-integrated B(xi, i?)-localized vortex-stretching 
terms transpires, 



(uj ■ V)it • <f>i uj dx ds = J —\uj(x,t)\ tpi(x) dx + 



\V(jj\ 2 d>j dx ds 



(4.7) 



-\uj\ 2 (((t)i) s + A(f)i) dx ds 
-\uj\ 2 (u ■ V4>i) dx ds, 
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for any t in (2T/3, T), and 1 < i < n. 

Denoting the time-averaged local vortex-stretching terms per unit mass associ- 
ated to the cover element B(xi, R) by VST Xi ^ Rt , 

(4.8) VST XitRtt = \J ^3 J(u-V)u-faujdxd s , 

the main quantity of interest is the ensemble average of {VST Xi]Ryt }™ =1 , 

n 

(4-9) (VST) R , t = - V 7ST x<At . 

n ^ — ' 

i=l 

Before stating the theorem, let us introduce the key macro-scale quantities, E , 
Po and <tq. Denote by Eqj time-averaged enstrophy per unit mass associated with 
the macro-scale domain B(0, 2R ) x (0, t), 

Eo,t — t [ 4? f \\u\ 2 <j>]! 2 dxds, 



tJo R 3 oJ 2' 

by Po.i a modified time-averaged palinstrophy per unit mass, 

If* If 1 1 f 1 

Po,t = -J -^3 y |Vw| 2 o ^rfs + --^3 y -|w(x,t)| ip (x)dx 

(the modification is due to the shape of the temporal cut-off r?), and by <7o,t a cor- 
responding Kraichnan-type scale, 



co,t 



Until now, there was no connection between the spatial macro-scale Rq and the 
global time scale T. At this point, it is convenient to assume i?o < \/T (in addition, 
without loss of generality, suppose that T < 1); in the case Rq > y/T, the proof can 
be modified similarly to the calculations in |DaG rll^2llDaGrl2-l| . 

Theorem 4.1. Let ube a global-in-time local Leray solution on R 3 x (0, oo), regular on 
(0, T). Suppose that, for some t e (2T/3, T), 

(4.10) C max{M| , i?| } a\ t < Rq 

where M = sup / \u\ 2 < oo, and C > 1 a suitable constant depending only on 

t Jb(0,2R ) 

the cover parameters. 
Then, 

(4.11) ^Po.t<(VST) R . t <CP .t 
for all R satisfying 

(4.12) C max{M(f , R^} cr| t < R < R . 

Remark 4.1. The macro-scale domain B(0, Rq) is placed at the origin for conve- 
nience only; it can be placed anywhere in K 3 . 
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Proof. Recall that 



(u> • V)it • fa w dx ds = J — \u)(x,t)\ ipi(x) dx + J j |Vw| <pi dx ds 



t , x 



o 



2 ,^;| 2 + dx ds 



(4.13) - j J^\uj\ 2 {u-Vfa) dxds, 

for any t in (2T/3, T), and 1 < i < n; the last two terms need to be estimated. 



For the first term, the properties of the spatiotemporal cut-off function 

2 rji > 2 p2 > 2 

3 - 3-^0 ^ 3 J 



setting pi = p2 = 3/4 - paired with the condition t > § T > |_R§ ^ |-R 2 yield 



(4.14) y y i|w| 2 ((^) s + Afa) dx ds < y y k| 2 0l /2 dx ds. 

The second term - the localized transport term - will be estimated similarly as 
in MGrZh061 ; the powers of the cut-off function fa will be distributed somewhat 
differently leading to a bit more precise estimate. 

Setting the cut-off parameters pi and p2 to 7/ 8, the following sequence of bounds 
transpires. 



^ ,oj| (it ■ Vfa) dx ds 



- C lll I (l w | 2 ^) 3/4 | u l (M 2 ^ 1/2 y ' >'■'■ d-s 



R Jo 

< C ^f\f \u\ A/ M 2 ^dxf" [J \u\ 2 <j>\' 2 dx)"* ds. 

(4.15) 

The first spatial integral is bounded as follows, 

\u\ 4/3 \u\ 2 fa dx 

1/3 



< (sup / |u| 2 dx) 2/3 ( / (|^|0 l 1/2 ) 6 dx) 

V » JB(xi,2R) ' V 7 

(4.16) < C (sup / |u| 2 dx) 2/3 ( / 1' 

(the last line by the Sobolev Embedding Theorem). 
Combining the bounds j4.15|l and | |4.16|| leads to 



dx 
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ri 



'1 



, \oj\ 2 (u ■ V 4>i) dx ds 
lo J 2 



<cifsup/ \u\ 2 dx) 1/2 (f [ \V(^ /2 Lj)\ 2 dxds) 3/4 ( f f ^(j)]' 2 dxds 

(4.17) 

< ^7 |v(#Ml' dxds + c ( SUP * J '«-ff> H2 *)' i, jf / m^**. 

Utilizing the commutator estimate (with p\ = p2 = 3/4) 
| \V(4u)\ 2 dx 



<2 / \Vcj\ 2 6 id x + C^ I \uj\ 2 6) l2 dx 



R 2 
(4.18) 

in the first term of the above inequality yields the final bound for the localized 
transport term, 



"i/ / \ Voj \ 2 (j>idxds + Crlp J J \lo\ 2 (/>1 /2 dxds 
Note that the factor 

R 

is scaling-invariant, and - consequently - the bound | |4.19l l is dimensionally cor- 
rect. However, for an arbitrary global-in-time local Leray solution it is not a priori 
bounded (it is a priori bounded, e.g., assuming a uniform-in-time bound on the 
L 3 -norm; this, however, automatically implies regularity [ISS03J). The best one 

can do in general is to simply write sup / \u\ 2 dx < Mq. 

s JB(xi,2H) 

Taking this into account, and applying the bounds ([4.141 and | |4.19| | in the ex- 
pression J4.13|l - describing the dynamics of the vortex-stretching term localized 
to the cover element B(xi, R) - leads to 



(4.20) 



(lo • V)it • 4>i u) dx ds = / — |w(:c,t)| ipi{x) dx + / / |Vw| 4>i dx ds + I L , 
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where 

\ T i\^lJ o J \Vu\ 2 <f>idxds + Cmax{MlR 2 }^ J J \u\ 2 ^ /2 dx ds. 

Ensemble-averaging l|4.20t and utilizing the inequality i4.6\ several times im- 
plies that as long as 

(4.21) C max{M| , i?f } o\ t < R , 

^Po,t<(VST) Rjt <CP 0jt 

for all R satisfying 

C max{M| , R% } o\ t < R<R - 

□ 

Remark 4.2. Suppose that T is the first (possible) singular time, and that the 
macro-scale domain contains some of the spatial singularities (at time T). This, 
paired with the assumption that u is a global-in-time local Leray solution implies 

lim cr ,t = 0; 

hence, the condition J4.10& in the theorem is automatically satisfied for any t near 
the singular time T. 
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